The Simplest Piston Problem I: Elastic Collisions 
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We study the dynamics of three elastic particles in a finite interval where two light particles are 
separated by a heavy "piston" . The piston undergoes surprisingly complex motion that is oscillatory 
at short time scales but seemingly chaotic at longer scales. The piston also makes long-duration 
excursions close to the ends of the interval that stem from the breakdown of energy equipartition. 
Many of these dynamical features can be understood by mapping the motion of three particles on 
the line onto the trajectory of an elastic billiard in a highly skewed tetrahedral region. We exploit 
this picture to construct a qualitative random walk argument that predicts a power-law tail, with 
exponent —3/2, for the distribution of time intervals between successive piston crossings of the 
interval midpoint. These predictions are verified by numerical simulations. 
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I. INTRODUCTION 



A classic thermodynamics problem is the adiabatic 
"piston" P], where a gas-filled container is divided into 
two compartments by a heavy but freely moving piston. 
The piston is clamped in a specified position and the 
gases in each compartment are prepared in distinct equi- 
librium states. The piston is then undamped and the 
composite system evolves to a global equilibrium. This 
simple scenario leads to surprisingly complex dynamics 
that are still incompletely understood, in both the cases 
where the two gases are elastic |^ H, ^ |^ and where 
they are inelastic 3- I^i the elastic system, the piston 
moves quickly to establish mechanical equilibrium where 
the pressures in each compartment are equal. Subse- 
quently, the piston develops oscillations that decayslowly 
as true thermal equilibrium is achieved 2, J, 0,1^. For 
the inelastic system, there is a spontaneous symmetry 
breaking in which the gas on one side of the piston gets 
compressed into a solid 0]. Surprisingly, this process is 
not monotonic, but rather, the piston undergoes oscilla- 
tory motion whose period grows exponentially with time. 

Given the complexities of these many-body problems, 
we instead investigate a much simpler version (Fig. 
a three-particle system ^ in the interval < a:; < 1 
consisting of two light particles of masses mi = — 1 
that are separated by a heavy piston of mass TO2 3> 1. 
All interparticle collisions and collisions between particles 
and the ends of the interval (henceforth termed walls) are 
elastic. We will develop a simple geometric approach and 
complementary numerical simulations to help understand 



the complex dynamical features of this idealized system. 
These results may ultimately be useful for understanding 
the many-body piston problem. 

An additional motivation to investigate the three- 
particle system is the connection to the collective behav- 
ior in one-dimensional (ID) fluids. The dimensional con- 
straint induces strong interparticle correlations that lead 
to anomalous transport properties. For example, heat 
conductivity is generally extremely large for ID fluids, 
while mass diffusion is exceedingly slow . An example 
of a ID fluid that exhibits such phenomenology is a gas 
of point particles with alternating masses 0. This fluid 
can be viewed a collection of three-particle subsystems, 
each similar to our idealized model. We therefore antici- 
pate that the dynamics of our three-particle system can 
shed light on anomalous collective phenomena that arise 
in ID fluids. 

In the next section, we outline the basic phenomenol- 
ogy of the three-particle system. Then in Sec. Ill, we 
map the trajectories of three particles on the line onto 
an equivalent elastic billiard particle that moves within 
a highly skewed tetrahedron, with the specular reflec- 
tion whenever the billiard hits the tetrahedron bound- 
aries [T3 . IT^ il4] . From this simple geometrical 
mapping, we can understand many of the unusual dy- 
namical properties of the system, as will be discussed 
in Sec IV. Perhaps the most unexpected feature is the 
long excursions of the piston close to the walls. By the 
billiard equivalence, we will argue, in Sec. V, that these 
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FIG. 1: The three-particle system — a piston and two light 
particles. 
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FIG. 2: (Color online) Trajectory of the piston for t < 400 for 
7712 = 100 and initial positions xi = 0.3, X2 = 0.5, and — 
0.7636; note the offset of the vertical axis. Inset: Trajectory 
for t < 10'', as well as the estimate X2{t) « Ei{t)/[Ei{t) + 
E3{t)] (see text). These data are averaged over a 400-point 
range for ease of visualization. 

long excursions can be understood as the collision point 
of the billiard in the tetrahedron undergoing a random 
walk. We will thereby find that the distribution of inter- 
val midpoint crossing times by the piston has a power-law 
tail with exponent —3/2. 

In an accompanying paper, we will consider the three- 
particle system when collisions between the light parti- 
cles and the walls are inelastic. Surprisingly, much of 
the phenomenology of this idealized three-particle sys- 
tem closely mirrors the complex dynamics that arises in 
the many-particle inelastic piston problem |^. 



II. PISTON MOTION 

To appreciate the basic phenomena, we show a typical 
piston trajectory obtained numerically in Fig. [3 for the 
case 1712 — 100. We focus on the case m2 — 100 because 
the overall simulation time scales as y^mj and there is 
already good scaling behavior for m2 = 100. However, 
we find qualitatively similar behavior for other values of 
7712 ^ 1- For concreteness and simplicity, we also con- 
sider the initial condition in which the light particles ap- 
proach a stationary piston at X2 = 1/2 with velocities 
vi = +1 and V3 = —1, so that the total momentum equals 
while the total energy E = 1. The initial particle po- 
sitions of the light particles are uniformly distributed in 
(0,1/2) and (1/2,1), respectively. Again, we also studied 
different initial speeds for the light particles and verified 
that our asymptotic results for the main quantities of 
interest do not qualitatively depend on the initial condi- 
tions. 

For the initial condition in which the light particles 
approach the piston with velocities ui = -t-1 and V3 = 
— 1, we find that after a short transient for t < 100, the 
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FIG. 3: (Color online) Particle energies as a function of time 
averaged over a 400-point range. 

piston settles into a quasiperiodic motion with period of 
T « 12. This time scale can be understood from simple 
arguments: if there is energy equipartition, the piston 
would have energy 1/3 and speed \v2\ = ^/2/3m^ « 0.08. 
Equipartition also implies that the typical spatial range 
of all three particles should be equal. These two features 
lead to a period of the order of 12, in agreement with the 
data. 

In the limit 7712 — > 00, it can be shown that the piston 
position obeys (see Ref. ^15] and the Appendix) 

(fx2{ts) _ Ai ^3 

corresponding to a particle moving in an effective po- 
tential well Vas{x) = ^[Aix^^ + yl3(l - x)'^]. Here 
ts = t/ y/m2 is a slow time variable that is a natural 
scale for the piston motion, and ^1^3 are initial condition- 
dependent constants. For total energy _E = 1, we nu- 
merically determine from Eq. ^ that the oscillatory 
period in slow time coordinates is T, « 1.285. Thus 
a piston with ?7i2 = 100 should oscillate with period 
T = Tsy/m2 ~ 12.85, in excellent agreement with sim- 
ulations (Fig. |2J|. Thus this effective potential picture, 
which formally applies in the limit 777,2 — 00, quantita- 
tively accounts for the short-time oscillations of a heavy 
but finite-mass piston. 

For t > 2000, however, a considerably slower and much 
less predictable large-amplitude modulation is superim- 
posed on the quasiperiodic oscillations (inset to Fig. [SI . 
When 7772 = 100, the piston eventually approaches to 
within of 0.05 of each wall. These long-time extreme ex- 
cursions are reflected in the time dependence of the par- 
ticle energies (Fig. EJ. For t > 2000, the piston energy 
fluctuates strongly and is phase locked with X2 (t) during 
the extreme excursions. Notice also that for t < 2000 the 
piston energy is consistently below its average long-time 
value, indicating the extent of the transient regime. 

We may alternatively estimate the piston position by 
mechanical equilibrium and basic thermodynamics. We 



write Piii (X Ei, where Pi and Ei are, respectively, the 
pressure and energy associated with particles i = 1,3, 
and £i is the length available to particle i. Assuming 
mechanical equilibrium. Pi = P3, and using ii = l — t^ — 
X2, we find X2(i) = Ei{t)/[Ei{t) + E3{t)]. This is very 
close to the numerical data for X2{t) (inset to Fig.|21); thus 
the piston excursions and the large energy fluctuations 
away from equipartition are closely connected. 

Finally, Eq. (Q) can be derived heuristically. The equa- 
tion of motion for the piston is m2V2 = F1—F3, where the 
overdot denotes the time derivative and Fi is the force ex- 
erted on the piston by particle i. Consider a time range 
large compared to the typical time between successive 
light particle bounces, but small compared to the time for 
the piston to move a unit distance. Then Fi « Api/ AU, 
where Api — —Avi is the momentum change of the piston 
after a collision with particle i, and Ati is the time be- 
tween successive bounces of the particle with the piston. 
When particle 1 with velocity vi , collides with the piston 
with velocity V2, the outgoing velocity of the former is 

1^2 

2v2 — vi in the limit 7712 ^ 1. Since V2 ^ 0{m2 ), 
we have Api « 2vi. Now Ati « 2ii/vi, where £1 — X2 
is the length of the subinterval that contains particle 1. 
Parallel results hold for collisions between the piston and 
particle 3. Thus 



•m2V2 = — 

X2 



1- X2 



(2) 



To obtain vi, note that after reflection from the left 
wall, particle 1 approaches the piston with velocity vi — 
2v2, so the net change in vi between successive collisions 
with the piston is —2v2. Thus the velocity of particle 1 
evolves according to vi ~ — 2v2/(2-^i/'yi) = —V1X2/X2, 
with solution vi cx l/x2- An analogous equation holds 
for V3. Using these results in Eq. Q gives Eq. 



III. BILLIARD MAPPING 

To help understand the unusual features of the par- 
ticle trajectories, it proves useful to follow conventional 
practice 10, 11, 12, 13, 14] and map the three-particle 
system into an equivalent effective billiard. To be gen- 
eral, suppose that the particles have masses mi, m2, and 
TO3, are located at < Xi{t) < X2{t) < 2:3 (t) < 1, and 
have velocities vi{t), V2{t), and V3{t). The trajectories 
of the three particles in the interval are then equivalent 
to the trajectory {xi{t), X2{t), X3{t)) of an effective bil- 
liard particle in the three-dimensional domain defined by 
the constraints < cci < 0:2 < CC3 < 1. For example, a 
collision between particle 1 and the left wall corresponds 
to the billiard ball hitting the boundary xi — 0, while 
a collision between particles 1 and 2 corresponds to the 
billiard hitting the boundary xi — X2, etc. 

Unfortunately, momentum conservation shows that 
collisions between the effective billiard and the bound- 
aries of the domain are not specular. Consequently, a 
naive analysis of successive billiard collisions becomes 




FIG. 4: The allowed tetrahedron (outlined by heavy lines) for 
the effective billiard particle in yt coordinates. The back and 
top planes of the tetrahedron are defined by j/i = and 7/3 = 
1, while the planes 1/2 — yiy/m2 and j/2 = yzy/rna correspond 
to 1-2 and 2-3 collisions. 



prohibitively cumbersome. However, a considerable sim- 
plification is achieved by introducing the "billiard" coor- 
dinates [13 111 111 111 111 



Hi = Xi y/rui 



w, = Vi ^ rrii 



1,2,3. (3) 



In these coordinates, the constraints xi < X2 and X2 < X3 
become 



m2 

yi < 2/2 : 

mi 



/TO3 

V2<V?. : 

1712 



while the constraints involving the walls are j/i > and 
2/3 < Vw3 = 1. As shown in Fig. ^ the allowed re- 
gion for the billiard is the interior of a highly skewed 
tetrahedron whose two acute interior angles are given by 
6 = tan^^ y^l/ rn2 ■ While this geometry may seem com- 
plicated at first sight, these coordinates ensure that all 
billiard collisions with domain boundaries are specular 
0, 0, Q , and this feature greatly simplifies the prob- 
lem. 

We now exploit this billiard mapping to characterize 
the motion of the piston in the original three-particle 
system. For a zero-momentum initial condition, the 
initial billiard trajectory lies within the shaded square 
2/2 = jJm2/2 (equivalent to X2 = 1/2 in the interval) in 
Fig. El If the first collision is between the piston and 
particle 1, the equivalent billiard first hits the front wall 
of the tetrahedron. Because of specularity, the billiard 
is refiected toward increasing 7/2 • Conversely, if the first 
collision is between the piston and particle 3, the billiard 
first hits the bottom wall and the reflected trajectory is 
toward decreasing j/2- 

The opposite effects of successive 1-2 and 2-3 colli- 
sions lead to the billiard persisting close to the shaded 
square y2 = y^TO2/2. However, once the billiard devel- 
ops a nonzero velocity in the 1/2 direction, the trajec- 
tory is unlikely to return to the initial square. Subse- 
quently, the billiard bounces back and forth primarily 
along the 1/2 direction in the tetrahedron, corresponding 
to the quasiperiodic oscillations in the interval shown in 
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Fig. 12 At still longer times, the billiard motion consists 
of unpredictable modulations that are superimposed on 
the quasiperiodic oscillations. The long-lived excursions 
of the piston near one end of the interval correspond to 
the billiard remaining close to one of the acute-angled 
ends of the tetrahedron in Fig. 01 

Another useful consequence of the billiard mapping is 
that we can also deduce in simple terms the probabil- 
ity distribution Tr2{x) for finding the piston at position 
a;2 = a; in the interval, or equivalently the probability for 
finding the billiard with coordinate j/2 — 'Xi^Jrrvi = z. 
This distribution can also be found by using the micro- 
canonical ensemble (see e.g., the first paper in Ref. 0)- 
If the billiard covers the tetrahedron equiprobably, then 
■K2{x) would be proportional to the area of the rectangle 
defined by the intersection of the plane y2 = z and the 
tetrahedron in Fig. ^ This mixing property is believed 
to occur in triangular billiards with irrational angles [T^ 
and also in various three-dimensional billiard geometries 
[TtI . Given that the angles of our tetrahedron generically 
are irrational except for particular values of TO2 , we ex- 
pect that billiard trajectories in this tetrahedron will also 
be mixing. 

From this mixing hypothesis, '!T2{z) is simply propor- 
tional to the area of the rectangle 2/2 — z in the tetrahe- 
dron. Now the length of the horizontal side of the rectan- 
gle is proportional to z/^m2, while the length of the ver- 
tical side is proportional to 1 — 2/ ^/rri2- Thus the rectan- 
gle area is proportional to z/^TO2(1 — zj ^Jrfvi) — 2:2(1 — 
X2)- Normalization of this probability fixes the propor- 
tionality constant and we thus obtain 'K2(x) — %x{\ — x) 
for the probability that the piston is located at x. Sim- 
ilarly, the position distribution of the light particles is 
found by computing the areas of the triangles defined by 
the intersection of the planes yi = x^fml (i = 1,3) with 
the tetrahedron. This leads to 'ni{x) — 3(1 — xY and 
"■3(2;) — 3a;^. 

We tested these predictions numerically and obtained 
excellent agreement between the above theoretical expec- 
tations and the simulation results. Notice that under the 
assumption of the billiard visiting all points in the tetra- 
hedron equiprobably, the probability of finding any par- 
ticle at a given position on the interval is a constant; that 
is, n(a;) = lY.i'^^i^) = 1- 
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FIG. 5: (Color online) The distribution of midpoint crossing 
time intervals P{5t) for ni2 — 2", with n = 4,6,8, and 10 
on a double logarithmic scale (lower to upper curves) for zero 
initial momentum of the system. For the case m = 64, we 
also show as open circles the distribution for random initial 
conditions (positions and velocities) of the light particles. For 
visibility, each successive curve is shifted vertically upward 
by 10"~*. Early-time oscillations do not appear for large 7712 
because the histogram bin is larger than interpeak spacing in 
Fig. ini Inset; The cutoff St as a function of m2. 




FIG. 6: (Color online) The short-time behavior of P{5t) for 
different 7712. 



IV. EXTREME EXCURSIONS 

To characterize the wanderings of the piston near the 
ends of the interval, we study the probability distribu- 
tion P{6t) to have a time interval St between successive 
midpoint crossings by the piston. A midpoint crossing 
corresponds to the equivalent billiard crossing the plane 
2/2 = ^7712/2. As shown in Fig. |31 P{St) decays as the 
power law {St)~^ over a significant time range. For the 
case of 7712 — 1024, the data for P{6t) versus St are 
quite linear on a double logarithmic scale for St in the 
range [1.8 x 10^, 3.1 x 10^]). We measure the slope to be 



fi = 1.5203 ± 0.0024. At longer times, the data have an 
exponential cutoff P{St) ^ e^''*/*'^, where St ^ with 
A = 2.14 ± 0.01 (inset to Fig. EJ. Correspondingly, the 

average time between crossings varies as {St) ^ rrS^ ■ 
From the relation X2{t) k. Ex{t) l\Ex{t) -I- (i)] de- 
rived in Sec.ni the piston crosses the midpoint whenever 
i?i(t) = Ez{t)\ thus P[St) can also be interpreted as the 
probability that the inequality Ex (t) ^ £'3 (t) persists for 
a time St. This long-time persistence of energy asym- 
metry is in agreement with previous simulations of IZ? 
binary fluids Q , in which light particles were reported to 
trap energy and release it very slowly. 
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FIG. 7: (Color online) Positions of the three particles for the 
system of Fig. |5| near t = 4000. Inset: Finer detail near 
t = 3997. 

The early-time sequence of peaks in P{St) is simply 
related to the half-period of the short-time piston os- 
cillations (and its resonances), = ^Tsy/rn2, where 
Tg « 1.285 is the slow-time period associated with a par- 
ticle in the effective potential of Eq. Thus the first 
peak of P{6t) should be at St w 2.6, 5.1, 10.3, and 20.6 
for 1712 = 16, 64, 256, and 1024, respectively, very close 
to the results in Fig. 

To understand the long-time power-law in P{St), we 
show in Fig. 13 the particle trajectories from Fig. [5] dur- 
ing the extreme excursion near t « 4000. This excursion 
is driven by a sequence of nearly periodic oscillations due 
to precisely orchestrated correlated motion of the lighter 
particles. Consider first the collisions between particle 
1 and the piston when the latter moves toward x = 0. 
There is a violent series of "rattling" collisions as the 
piston first approaches x — and ultimately is reflected 
|14L . In the limit m2 ^ 1 , these rattling collisions are 
equivalent to the piston having a nearly elastic reflection 
from the wall. After this rattling collision, the piston is 
met by particle 3 whose momentum is of a similar mag- 
nitude, but opposite to that of the piston. Thus after 
a few collisions between the piston and particle 3 (seven 
such collisions in Fig. [7^ the piston is reflected back to- 
ward X = 0, where the rattling between particle 1 and 
the piston recurs. 

V. EFFECTIVE RANDOM WALK 
DESCRIPTION 

To determine P{6t) from this descriptive account of 
the rattling collisions, we consider a reduced problem in 
which the fastest degrees of freedom associated with par- 
ticle 1 are integrated out. As we shall show, the slower 
degree of freedom associated with the piston can then 
be described qualitatively by an effective random walk. 
This then allows us to deduce the statistics of the time 



between successive piston crossings of the interval mid- 
point. 

For the piston to persist near the left wall, the col- 
lisions between the piston and particle 3 must be close 
to periodic. A deviation from periodicity occurs because 
the net effect of the rattling between particle 1 and the 
piston is a slightly inelastic collision. We now estimate 
the departure from elasticity in these rattling collisions 
and we then use this result to estimate the duration of 
the resonance between the piston and particle 3. 

In billiard coordinates, the rattling collision can be rep- 
resented as the effective billiard entering a narrow wedge 
of opening angle 9 = tan~^(l/^/m2) [Fig. Ela)] that is 
the projection of the tetrahedron onto the yi-y2 plane. 
Because each collision of the billiard with the wedge is 
specular, the ensuing rattling sequence is equivalent to a 
straight trajectory in the periodic extension of the wedge 
[Fig.lHfa)]. Each collision is alternately particle-particle 
and particle- wall, so that the identity of periodically ex- 
tended barriers alternates between pp and pw. The rat- 
tling sequence ends when the billiard trajectory no longer 
crosses a wedge boundary. The crucial point is that the 
flnal billiard velocity vector deviates by no more than an 
angle with respect to the two rays that deflne the last 
wedge. 

Suppose that the initial velocity vector is tT*^'-* = 
{vi,V2) = (0,-1), corresponding to w^'-* = {wi,W2) — 
(0, —^7712). If the final billiard trajectory is parallel to 
a pw boundary in Fig. (Hfa), then v^^^^ is (0,-1-1). This 
corresponds to a rattling sequence in which particle 1 
begins and ends at rest and the piston is elastically re- 
flected. Conversely, if the flnal trajectory is parallel to a 
pp boundary, then 

Wf = ^m2/(l + 1712) (-1, -\/m^) 

(note that wf — w'j). Translating to original coordi- 
nates, the minimum flnal piston speed is 1112/ (1 -t- m2). 
Therefore rattling collisions lead to a flnal piston veloc- 
ity that lies within the narrow range (1 — l/2m2, 1) for 
TO2 > 1- 

From this deviation from elasticity, we determine the 
time needed to disrupt the resonance between the pis- 
ton and particle 3. Consider the two-particle system 
consisting of the piston and particle 3 with initial ve- 
locities (w2,i'3) — (l/™2,— 1) and with X2 = X3 = x = 
1/(1-1- 1712)- This resonant starting state ensures that 
the two particles hit the opposite ends of the interval si- 
multaneously and then meet again at x = 1/(1-1- 7112) 
with {v2,V3) = (1/to2,— 1) when all collisions are elas- 
tic. In billiard coordinates this periodic motion trans- 
lates to a singular trajectory in the projection of the 
tetrahedron onto the 2/2-2/3 plane [Fig. [SJb)]. A 2-3 col- 
lision corresponds to the billiard hitting the hypotenuse 
of the resulting triangle perpendicularly and the simul- 
taneous collision of the particles with the interval ends 
corresponds to the billiard hitting at the right-angle cor- 
ner of the triangle. If the initial position slightly deviates 
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FIG. 8: (a) The wedge that represents coUisions between par- 
ticle 1, the piston, and the left wall in billiard coordinates 
(thick lines). Lighter lines show the periodic extension of 
wedge. The dashed straight line is periodic extension of the 
billiard trajectory, (b) The billiard after projection from the 
tetrahedron onto the j/2-1/3 triangle. Shown is the periodic 
trajectory when the piston and particle 3 have equal and op- 
posite momenta and meet at x — 1/(1-1- 1712)- A two-cycle 
arises (dashed) when the collision point deviates within the 
thick segment of the hypotenuse while the particle momenta 
remain equal and opposite. 



from X = 1/(1-1- TO2) while the initial velocities are on 
resonance, then the subsequent motion is simply a two- 
cycle, as indicated in the figure. 

Now consider the influence of particle 1 on this reso- 
nance. Due to the slight inelasticity of the effective col- 
lision between the piston and the left wall, the speed of 
the reflected piston changes stochastically by the order 
of m^^. Consequently, the 2-3 collision point shifts from 
X = 1/(1 + TO2) to a: = 1/(1 + ma) ± 0{l/ml). In bil- 
liard coordinates, the return trajectory in Fig. |Sfb) is 
not exactly parallel to the initial trajectory and the colli- 
sion point on the hypotenuse moves stochastically by an 

3/2 

amount of the order of . When this collision point 

moves outside the thick line in Fig. |S1 the resonance be- 
tween the piston and particle 3 terminates and the piston 
crosses the interval midpoint shortly thereafter. 

Thus the collision point on the hypotenuse undergoes a 

— 1/2 

random walk on an interval of length 0(m2 ) with one 
end absorbing (open circle in Fig. (S)) and the other end 
reflecting (solid dot). The probability that the billiard 
remains in this interval up to time 6t therefore scales 
as P{St) ~ St~^^^ until an exponential cutoff because 
of the finiteness of the interval |0. The cutoff time 

should be L^/D, where L ~ 0(m^^^^) is the interval 

length and D cx [0(m,2^^^)]^ is the diffusion coefficient 
associated with individual random-walk steps of length 
7712 '^^^ ■ This leads to a cutoff time St m^, consistent 
with the simulation result St ~ rn^'^^ shown in Fig. |S1 



VI. DISCUSSION 

We introduced a toy version of the classic piston prob- 
lem in which a massive particle (the piston) separates a 
finite interval into two compartments, each containing a 
single light particle. In spite of its simplicity, the dy- 
namics of this three-particle system is surprisingly rich. 
When all collisions are elastic, the piston undergoes com- 
plex motion, with short-time quasiperiodic behavior and 
seemingly chaotic behavior at long times. The early-time 
can be understood in terms of the piston oscillating in 
an effective potential well Vcff — [Aix^"^ -f ^3(1 — a;)~^]. 
To understand the long-time behavior, we mapped the 
motion of three particles in the interval onto that of an 
effective equivalent billiard in a tetrahedral domain. We 
then used geometric methods that help explain some of 
the anomalous dynamical features of the piston. 

At long times, the piston moves in an apparently un- 
predictable fashion, with long-lived excursions close to 
the ends of the interval during which large departures 
from energy equipartition occur. We quantified these 
extreme excursions by studying the distribution of time 
intervals for the piston to cross the interval midpoint. 
This distribution has a power-law decay over a wide time 
range, with exponent —3/2. We argued that this phe- 
nomenon can be recast as a first-passage problem of a 
random walk within a finite interval, leading naturally 
to the above exponent value. 

Although individual trajectories of the piston seem un- 
predictable, average properties are considerably simpler. 
We found a simple form for the probability distribu- 
tion 'Ki{x) for finding particle i at position a;; namely, 
7ri(x) = 3(1 — x)^, 7r2(a;) = 6x(l — x), and 7r3(a;) = Zx^. 
These forms are a direct consequence of the billiard tra- 
jectory being mixing in the tetrahedron. 

The three-particle system studied in this work is clearly 
oversimplified to faithfully model a many-particle piston 
system in three dimensions. Nevertheless, the methods 
developed here may prove useful in understanding few- 
particle elastic or granular systems and may help suggest 
new approaches to deal with many-particle systems in 
higher dimensions. 
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APPENDIX: EFFECTIVE POTENTIAL FOR AN 
INFINITE-MASS PISTON 



Here we show that Eq. |^ governs the piston motion 
in the Hmit m2 — > c» by speciahzing the general result of 
Sinai jl^ to the three-particle elastic system in the inter- 
val. We assume a initially at X2(0) = 1/2 with ^2(0) = 0, 
and unit-mass particles 1 and 3 starting at random posi- 
tions to the left and right of the piston, respectively, with 
velocities ui(0) = -f 1 and ^3(0) = —1. Energy conserva- 
tion implies that |w2(i)| < ■\/2/m2. It is then natural to 
define a slow time variable, ts — tj ^mi-, such that the 
piston velocity is 0(1) in this time scale. 

Consider an infinitesimal slow time interval [t^ , -I- (5] 
during which the piston moves a distance 0((5), while the 
number of 1-2 and 2-3 collisions is 0(^7712). Let fc index 
each piston collision; we define this collision index to run 
from A;_ -|- 1 to A;+ in [t^ , -\- S\. The total number of 
collisions experienced by the piston in this time range is 
N — fc+ — fc_. The particle velocities just before each 
collision with the piston are given by 



V2(fc) = (l-e)w2(fc-l) + et;,(fc-l), (A.l) 
v^(k) = (e-l)w,(fc-l) + af2(fc-l), (A.2) 



where e = 2/(1 -I- TO2), = e'7i2, i = 1,3, and w(fc) is 
a particle velocity just before the {k -I- l)st piston col- 
lision. For large piston mass recoUisions do not occur, 
that is, light particles always hit a boundary before col- 
liding again with the piston. Therefore, ui(fc) > and 

v^ik) <o\f ke [k_,k+-i]. 

Next we iterate the first term in Eq. (|A.ip to write 
(fc+) = ^2(^3 + S) in terms of W2(fc-) = V2{ts). Let ni2 
and 7123 be the number of 1-2 and 2-3 collisions in the 
sequence A:_ -I- 1, . . . , fc_|-, respectively, with N — ni2 + 
77.23. For i e [1,7112], we define ci{i) = fc if and only if 
the 7th 1-2 collision corresponds to the fcth collision in 
k- + 1,...,A;+, so that ci{i) £ [1,-/V] and similarly for 
csO), with j e [1,7723]. With these definitions, HA.1|) 
gives 



V2{k+) = {l-efv2{k^) 



(A.3) 



e^(l-e)^-=i« 7'i[fc_+ci(i)-l] 

i=l 

"23 

6^(l-e)^-^^W 7;3[A:_-KC3(7)-1]. 



The piston velocity in the slow time variable is W2{ts) = 
^ = ^Jrn2 ^2 (i s ) ■ To derive a closed equation for 7^2 {tg ) , 
we first take the limit 7772 — > 00 and then (5 — s- 0. Using 



the definition of 702 (is) and Eq. ljA.3|) . we find 

W2{ts + 5) - W2{ts) = [(l-e)^-l]7i;2(i.) 
+ey/T^[ J2 [(1 - e)''-^'^'^ - 1] vi[k. + ci(7) - 1] 



. g [(1 _ ,f-csU) _ 1] „3[fc_ + C3(j) - 1] 



"12 



- ^ V, [k_ + Ci (7) - 1] + £ 7^3 [fc- + C3 U) - 1] } 



(A.4) 



We expand this expression for 7772 00, taking into ac- 
count that e ~ 0{m2^) and 771^3 ~ 0{y^m2), to obtain 



"12 

W2{ts + S)-W2its) = ey/rn^^ ^ vi[k^ + ci{i) - 1] 

i=l 

"23 / -, 

^7;3[fc-+C3(jVl][ + — 



3=1 



/7772 



(A.5) 



Because the large piston mass causes the light particle 
velocities to change only slightly in the slow time interval 
[tg , is -f (5] , we can write 



Wi[fc- + Ci(fc) - 1] ~ riiViits 



(A.6) 



fc=i 



for i — 1,3, with Vi(ts) = Vi{k-), and where correction 
terms vanish as (5 ^ Within this approximation 

of nearly constant light-particle velocities, the unsealed 
time intervals between successive 1-2 and 2-3 collisions 
are 2x2{ts)vi'^{ts) and —2[\—X2{ts)]v^'^{ts)^ respectively. 
Thus 



ni2 



2x2{ts) 



n23 



2{l-X2{ts)) ' 



with vi{ts) > while 7;3(is) < 0. Using these results 
Eq. ETS)) . we obtain, in the asymptotic limit, 



X2{ts) [1-2:2(^5 



(A.7) 



We now derive the equation of motion for Vi{tg), i = 
1,3. Here we consider only particle 1, since the derivation 
for particle 3 is analogous. Let us introduce a new index 
q G [1,7712], such that vi{q) and V2{q) are the velocities 
of particle 1 and the piston just before the {q + l)st 1- 
2 collision (notice a subtle difference from the previous 
notation; between the gth and the (g-l-l)st 1-2 collisions, 
the piston may collide one or more times with particle 
3) . From Eq. HA.2|) we have 

7;i(g + 1) = (1 - e) vi{q) - a V2{q) ■ 
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Notice the extra minus sign in this equation compared to 
Eq. (|A.2|I to account for the reflection of the hght particle 
off the wall. Iterating this equation, we find 

" 12- 1 

viin,2) = (1 - 6)"- v,iO) - a ^ (1 - e)"--9-i ^^(g), 

9=0 

where now Wi(ni2) = vi{ts + 6) and wi(0) = vi{ts). Tak- 
ing the limit m2 — > oo now yields 



7112 — 1 



{ts +S)~ vi{ts) = -2 ^ V2{q) + O 



9=0 



1 



/m2 



We now write V2{q) as ^2(0) + J2l=i [^2(fc) ~ V2{k ~ 1)] . 
Therefore, 



"12-1 9 



^ V2{q) ^ ni2V2{0) + J2 I][«2(fc)-^;2(fc-l)]. 

9=0 9=0 fc=l 

(A.f 



The double sum in Eq. ljA.8|) can be demonstrated to be 
0{6^) [l^, so it is negligible in the 6^0 limit. Therefore 



Vi{ts +S) - Vi(ts) 



X2{ts) 



V2{ts) + 0{8'') . 



Finally, for TO2 00 and (5 — > we find 



dvi {ts 
dts 

dts 



Vljts) dx2 
X2{ts) dts 

vsjts) dX2 

1 - X2{ts) dts 



(A.9) 



Integrating these equations yields Vi{ts) = Bi/x2{ts) and 
VTiits) = B^/^l — X2{ts)], with Bi^3 constants which de- 
pend on the initial condition. Using these solutions in 
Eq. ljA.7|) we finally arrive to the piston equation of mo- 
tion given in Eq. l(T]l. 
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